Europhysics Letters 



PREPRINT 



Integrals over Products of Distributions 
from Manifest Coordinate Invariance of 
Perturbation Expansions of Path Integrals in Curved Space 

H. Kleinert(*)^ and a. Chervyakov(**) ^ 
Freie Universitdt Berlin 
Institut fiir Theoretische Physik 
Arnimallee 14, D-14-195 Berlin 

PACS. 3.65.-W - Quantum mechanics. 

PACS. 2.90.+P - Other topics in mathematical methods in physics. 
PACS. 2.30.Qy - Integral transforms and operational calculus. 



Abstract. - We show that the requirement of manifest coordinate invariance of perturbatively 
defined quantum-mechanical path integrals in curved space leads to an extension of the theory of 
distributions by specifying unique rules for integrating products of distributions. The rules are 
derived using equations of motion and partial integration, while keeping track of certain minimal 
features stemming from the unique definition of all singular integrals in 1 — e dimensions. 
Our rules guarantee complete agreement with much more cumbersome calculations in 1 — e 
dimensions where the limit e — > is taken at the end. In contrast to our previous papers where 
we solved the same problem for an infinite time interval or zero temperature, we consider here 
the more involved case of finite-time (or non-zero temperature) amplitudes. 



Introduction. - Until recently, a coordinate-independent definition of quantum mechan- 
ical path integrals in curved space existed only in the time-sliced formulation This is in 
contrast to field-theoretic path integrals between two and four spacetime dimensions which 
are well-defined in continuous spacetimes by perturbation expansions. Initial difficulties in 
guaranteeing coordinate independence were solved by 't Hooft and Veltman ||^ using dimen- 
sional regularization with minimal subtractions (for a detailed description of this method see 
the textbook Q). Coordinate independence emerges after calculating all Fcynman integrals 
in an arbitrary number of dimensions d, and continuing the results to the d = A. Infinities 
occuring in the limit are absorbed into parameters of the action. 

In contrast, and surprisingly, numerous attempts |]7|-]l5|| to define the simpler quantum 
mechanical path integrals in curved space by perturbation expansions encountered problems 
in evaluating the Feynman integrals. Although all physical properties are finite and uniquely 
determined by Schrodinger theory, the Feynman integrals in the expansions are highly sin- 
gular and mathematically undefined. When evaluated in momentum space, they yield dif- 
ferent results depending on the order of integration. Various definitions chosen by earlier 
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authors |[^-[l5|] were not manifestly coordinate- independent, and this could only be cured by 
adding coordinate-dependent "correction terms" to the classical action — a highly unsatisfac- 
tory procedure violating the basic Feynman rule that physical amplitudes should consist of 
a sum over all paths of phase factors e*"^ whose exponents contains only the coordinate- 
independent classical action along the paths. 

The first satisfactory perturbative definition of path integrals in curved space was found 
only recently by us [^-Q. The results enabled us to set up simple rules for treating integrals 
over products of distributions in one dimension to ensure coordinate invariance |4). These 
rules were derived for path integrals on an infinite time interval or zero temperature, where 
we could apply most directly the dimensionally continued integration rules of 't Hooft and 
Veltman |^ in momentum space. 

In a recent pape r j|l6|] , the authors of and (l4| have adapted the methods developed in 
our first two papers to the calculation of finite-time amplitudes (see also [0|l8| ) . In doing 
this they have, however, not taken advantage of the great simplifications brought about by 
the developments in our third paper [Q which make explicit evaluations of Feynman integrals 
in d = 1 — e dimensions superfluous. The purpose of the present work is to show how this 
happens. Thus we shall derive rules for calculating integrals over products of distributions 
which automatically guarantee coordinate independence. All integrals will be evaluated in 
one dimension, after having been brought to a regular form by some trivial manipulations 
which require only a small residual information on the initial 1 — e -dimensional nature of the 
Feynman integrals. 

Consider the short-time amplitude of a particle in curvilinear coordinates 

{qbn\qaTa)= I Pg(r)V^e--^[«l, (1) 



where A[q] is the euclidean action of the form 

'1 



A[q] = / dr 



^g,,{q{T))q^T)q^T) + V{q{T)) 



(2) 



The dots denote r-derivatives, gij{q) is a metric, and g = detg its determinant. The path 
integral may formally be defined perturbatively as follows: The metric gij{q) and the potential 
V{q) are expanded around some point near qa and qb in powers of Sq^ = q^ — q^. After 
this, the action A[q\ is separated into a free part Avi[qo\5q\ = (1/2) j^^ dr gij{qo)q^q^ , and an 
interacting part Antfeo; 5q] = A[q] - Ao[qQ;dq]. 

A simply curable ultraviolet (UV) divergence problem is encountered in the square root 
in the measure of functional integration in (|l|). Taking it into the exponent and expanding in 
powers of Sq, it corresponds to an effective action 

^Va = -^(0) I ' dTlog[g{qo + Sq)/g{qo)], (3) 

which contains the (5-function at the origin (5(0). This infinite quantity represents formally the 
inverse infinitesimal lattice spacing on the time axis, and is equal to the momentum integral 
5{0) = Jdp/{2'K). With (|), the path integral ^ takes the form 

(qbUlqaTa) = / I?g(r)e--^M--^v^M = / p(5g(r) e--^[*+*'?l--^v^['f'+''«l . (4) 



The main problem arises in the expansion of the amplitude in powers of the interaction. For 
simplicity, we shall set = 0, Tb = /?, as in thermodynamics and assume qb = qa = qo = 0. 
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Performing all Wick contractions, the origin to origin amplitude (0/3|0 0) is expressed as a 
sum of loop diagrams. There are interaction terms involving Sq^Sq"' which lead to Feynman 
integrals over products of distributions. The diagrams contain four types of lines representing 
the correlation functions 

A{t,t') ^ {6qiT)6q{T')) ^ , A-(r,r')EE(<5g(r)5Q(r'))= , 

■AiT,T') ^ {5q{T)6q{T')) = , ■A{T,T')^{dq{T)Sq{T'))= (5) 

The right-hand sides show the line symbols to be used in Feynman diagrams. 

The first correlation function A(r, r') = A(r', r) is determined by the free part ^ol^o; ^o] 
of the action. It is the Green function of the equation of motion 

■A(t,t') = A-(T,r') = -5(r-r'), (6) 

satisfying the Dirichlet boundary conditions 

A(0, r') = AiP, t') = 0, A(r, 0) = A(t, P) = 0. (7) 

Explicitly, it reads 

A(r, r') = A(r', r) = i [-e{r - t'){t - r') + t + t'] - (8) 

where e(r — r') is the antisymmetric distribution which is equal to ±1 for t ^t' and vanishes 
at the origin. 

The second and third correlation functions A (t, t') and A(t, t') are 

■A (r, r') = -ie(r -r') + ^^j, A(r, r') = ^eir - r') + ^ ~ ^ = A (r', r) , (9) 

with a discontinuity at r = r'. Here and in the following, dots on the right and left of A(t, t') 
denote time derivatives with respect to r and r', respectively. 
The fourth correlation function A'(t, r') is simply 

■A{t,t')^6{t-t')-1/P. (10) 

The 5-function arises from the derivative 6{t — r') = e(r — t')/2. Its value at the origin must 
be equal to the prefactor ^(0) of the effective action of the measure to cancel all ultraviolet 
(UV) infinities. Note the close similarity of ( p^ ) to the equation of motion (^). 

The difficulty in calculating the loop integrals over products of such distributions is best 
illustrated by observing the lack of reparametrization invariance of the path integral of a 
free particle in n-dimensional curvilinear coordinates first done by Gervais and Jevicki p|, 
Salomonson and recently also by Bastianelli, van Nieuwenhuizen, and collaborators |p|-p^|. 
The basic ambiguous integrals causing problems arise from the two-loop diagrams 



O : hi^ I I dTdT'A(T,r') A(r,r')A-(r,r'), (11) 
Jo Jo 

/i5= r /^dTdT'A(r,r')A-2(T,T'). (12) 



It is shown in Appendix A that the requirement of coordinate independence implies that these 
integrals have the values 

hi = /3/24, /fs = -/3/8, (13) 
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where the superscript R denotes the finite part of an integral. 

Let us demonstrate that these values are incompatible with partial integration and the 
equation of motion (^. In the integral (11), we use the symmetry 'A (t, r') = A'(r, r'), apply 
partial integration twice taking care of nonzero boundary terms, and obtain on the one hand 



Ji4 - - / dTdT'-A{T,T')—[A^{T,T')]^-- drdr' A^t.t')-A{t,t') 

^ Jo JO ^ Jo JO 



On the other hand, we apply Eq. ( [lO|) and perform two regular integrals, reducing /14 to a 
form containing an undefined integral over a product of distributions: 



/i4 = drdr' ■A{t,t')A{t,t')5{t -t')-\ f dr dr' 'A {t,t') A {t,t' 

JO JO P JO JO 



/3 rP 



P rP 



Jo 



drdr' 



-le'ir-r')6ir-r') 



/3 



(15) 



A third, mixed way of evaluating /14 employs one partial integration as in the first line of 
Eq. (|l^), then the equation of motion (||) to reduce /14 to yet another form 



/l4 = X / / dTdT'A'{T,T')d{T-T')^ 

^ Jo JO 

I rP rP 1 rP 

= - / dTdT'e^iT~T')S{T-T') + - dTA^{T,T) 

° Jo Jo 2 Jo 



f3 



We now see that if we set 



drHr)r5{T) 



(16) 



(17) 



the last two results ( |l6| ) and ( [T^ ) coincide with the first in Eq. (|l4[). The definition ([l7| ) is 
obviously consistent with partial integration if we insert S{t) = e(r)/2: 

J dr[eiT)r 6ir) = \ j dr[.(r)]2 ^(r) = \ j dT^[e{T)f = 1 (18) 

While the integration rule ( [l7| ) is consistent with partial integration and equation of motion, it 
is incompatible with the requirement of coordinate independence. This can be seen from the 
discrepancy between the resulting value /14 = /3/12 and the necessary ([T^). This discrepancy 
was compensated in Refs. [f7|~|l5| by adding the above-mentioned noncovariant term to the 
classical action. 

A similar problem appears with the other Feynman integral (^2|). Applying first Eq. 
we obtain 



j-P rP 2 /■'^ 1 /■'^ f'^ 

/i5- dTdT'A{T,T')S^{T^T')~-J^ dTA{T,T) + ^j^ drdr' A{t, t' ) . 



P rP 



(19) 
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For the integral containing the square of the (5-function we must postulate the integration rule 



/ 



dT[<5(r)]V(r)^5(0)/(0) (20) 



to obtain a divergent term 

ltt = SiO)jyrA{r,r)^S{0)^. (21) 

proportional to (5(0) compensating a similar term from the measure. The remaining integrals 
in (19) are finite and yield the regular part of J15 



/il = -f • (22) 



In another calculation of /15, we first add and subtract the UV divergent term, writing 

/l5 = r r dTdT'AiT,T') [■A^iT,T')-6^iT-T')] +6(0)^. (23) 

Jo Jo " 

Replacing S^{t — t') by the square of the left-hand side of the equation of motion (^, and 
integrating the terms in brackets by parts, we obtain 

r/3 rP 



/« = / / drdr'A(T,T')[-A-2(r,r')- A-2(r,r')] 

Jo Jo 

= drdT'[-A(T,T')A(T,T')A(T,T')- A(r,T')A-(r,r')'A-(r,r')] 

Jo Jo 

- drdr' [~A^{t,t')A-{t,t')- A{t,t')A{t,t')A-{t,t')] 

Jo Jo 

= -/l4+ / dTdT'A^{T,T')A-{T,T') = -Iu-f3/6. (24) 

Jo Jo 



The value of the last integral follows from partial integration. 

For a third evaluation of J15 we insert the equation of motion (^ and bring the last integral 
in the fourth line of (p3) to 



f drdr' A'^{t,t')5{t-t') = -I3 - [ dr e'^{T)S{T) + — 
Jo Jo [4 J 12, 



(25) 



All three ways of calculation lead to the same resuh /^l = -/3/4 usmg the rule (0). This, 
however, is again in disagreement with the coordinate- invariant value in Eq. (|l3|). Note that 
both integrals In and //I are too large by a factor 2 with respect to the necessary ( |l3| ) for 
coordinate invariance. 

How can we save coordinate invariance while maintaining the equation of motion and 
partial integration? The direction in which the answer lies is suggested by the last line of 
Eq. ([l6|): we must find a consistent way to have an integral 

dT[e{T)fd{T)^0, (26) 
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instead of ([l^), which means that we need a reason for forbidding the apphcation of partial 
integration to this singular integral. For the calculation at the infinite time interval, this 
problem was solved in our previous papers with the help of the dimensional regulariza- 
tion, carried to higher orders in Refs. The extension of our rules to the short-time 

amplitude considered here is straightforward. It can be done without performing any of the 
cumbersome calculations in 1 — e-dimension. We must only keep track of the essential features 
of the structure of the Feynman integrals in arbitrary dimensions. For this we continue the 
imaginary time coordinate t to a rf-dimensional spacetime vector t ^ = {t, , . . . , x'^~^). 
In d = 1 — e— dimensions, the correlation function reads 

A(r,x; r',x') = J ^e^H^-^'^ AM,r'). (27) 

Here the extra e-dimensional space coordinates x are assumed to live on infinite axes with 
translational invariance along all directions. Only the r-coordinate lies in a finite interval 
< T < /3, with Dirichlet boundary conditions for (|27|). The one-dimensional correlation 
function A^{t,t') in the integrand has a mass uj = |k|. It is the Green function on the finite 
r-interval 

- •A^(r, r') + Lu^AUr, t') = 5{t - r'), (28) 
satisfying the Dirichlet boundary conditions 



Explicitly, it reads Q] 



A^(0,t) = A^(/3,r) =0. (29) 

A f-. sinh^(/3- r>) sinhtJT< 
Auj{t, t ) = — — ^ , (30) 

where r> and r< denote the larger and smaller of the imaginary times r and t', respectively. 
In d dimensions, the equation of motion (H) becomes a scalar field equation of the Klein- 



Gordon type. Using Eq. (28), we obtain 



^^A(t,x; t',x') = A^^(r,x; r',x') =-A(t,x; t',x')+ xxA(t, x; t', x') 



i!^,.k(x-x') [•A^(r,r')-c.2A^(r,r')] 



= -5{T-T')5'^'\-K-y:!) = -5'^'^\x~x'). (31) 

The important observation is now that for d spacetime dimensions, perturbation expansion 
of the path integral yields for the second correlation function A' (r, r') in Eqs. ( [Tl| ) and (|2|) the 
extension ^Ai,{x,x' ). This function differs from the contracted function ^A^(a;,a;'), and from 
^^A(a;,a;') which satisfies the field equation (|3l|). In fact, all correlation functions A'(r, r') 
encountered in the diagrammatic expansion which have different time arguments always have 
the c?-dimensional extension ^Au(x,x'). An important exception is the correlation functions 
at equal times A(t, t) whose d-dimensional extension is always^ A^(a;, x), which satisfies the 
equation ( p^ in the e O-limit. Indeed, it follows from Eq. ( p7| ) that 

^,A^{x,x) = j [A^{t,t)+u:^A^{t,t)\ . (32) 

With the help of Eq. (^) , the integrand in Eq. (|3^) can be brought to 
.». / X 9. / N t-^^x UJ cosh uj{2t — 13) 

A.(r,r) +.^A.(.,r) = ^(0) . (33) 
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Substituting this into Eq. (|32D, we obtain 

^A,,{x,x)^S^''\x,x)~r . (34) 
The integral P is calculated as follows 

#fc wcoshw(2r-/3) I , ^ cosh 2(1 - 2r//3) 

— dzz — (35) 



{2ttY sinhw/3 (3 (27r/3)^ Jq sinh 

^(^+')-[C(. + l,l-;) + C(e + l/^ 



where 5'^ = 27r^/^/r(£/2) is the surface of unit sphere in e dimension, and T{z) and ({z, q) are 
gamma and zeta functions, respectively. For small e — s- 0, they have the limits t^{e + l,q) 
1/e — V'('z), and r(e/2) 2/e, so that ^ 1/(3, proving that the d-dimensional equation (Q) 
at coinciding arguments reduces indeed to the one-dimensional equation (0). The explicit 
d-dimensional form will never be needed, since we can always treat ^tA^(x,a;) as purely one- 
dimensional objects 'A'(r, r), which can in turn be replaced everywhere by the right-hand side 
6{0)-l/p oify. 

We now show that by carefully keeping track of the different contractions of the deriva- 
tives, we obtain a consistent calculation scheme which yields results equivalent to assuming 
integration rule ( p6|) in the calculation of /14 and /15, thus ensuring coordinate independence. 
The integral (^l]) for /14 is extended to 

Im= f f d''xd''x\A{x,x')A,{x,x')^A,{x,x'), (36) 



and the different derivatives on ^A^{x,x') prevent us from applying the field equation (|3l|), 
in contrast to the one-dimensional calculation. We can, however, apply partial integration as 
in the first line of Eq. (fih, and arrive at 



Ifi ^ -IJ J d''xd^'x'Al{x,x')A^^{x,x'). (37) 

In contrast to the one-dimensional expression (p[), a further partial integration is impossible. 
Instead, we apply the field equation (^l]), go back to one dimension, and apply the integration 
rule (|2^) as in Eq. ( [l^ ) to obtain the correct result /14 ~ /3/24 guaranteeing coordinate 
invariance. 

The Feynman integral (p^ for /15 is treated likewise. Its d-dimensional extension is 



4 - y y d''xd''x'A{x,x') lA,{x,x')Y . (38) 

The different derivatives on ^A^ (x, make it impossible to apply a dimcnsionally extended 
version of equation (|l^ as in Eq. (|9|). We can, however, extract the UV divergence as in 
Eq. ( |2^ ) , and perform a partial integration on the finite part which brings it to a dimensionally 
extended version of Eq. (|^ : 

I?5 = -hi + J d''xd''x'Alix,x')A^^{x,x'). (39) 

On the right-hand side we use the field equation (pi]), as in Eq. ([25|), return to d = 1, and use 
the rule ( [2^ ) to obtain the result //| = — /14 — /9/12 = — /3/8, again guaranteeing coordinate 
independence. 
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Thus, by keeping only track of a few essential properties of the theory in d dimensions 
we indeed obtain a simple consistent procedure for calculating singular Feynman integrals. 
All results obtained in this way ensure coordinate independence. They agree with what we 
would obtain using the one-dimensional integration rule ( p6| ) for the product of two e- and 
one (5-distribution. 

Our procedure gives us unique rules telling us where we are allowed to apply partial inte- 
gration and the equation of motion in one-dimensional expressions. Ultimately, all integrals 
are brought to a regular form, which can be continued back to one time dimension for a direct 
evaluation. This procedure is obviously much simpler than the previous explicit calculations 
in d dimensions with the limit d 1 taken at the end. 

We now apply this procedure to the perturbation expansion of the short-time amplitude 
of a free particle in curvilinear coordinates. 



Perturbation Expansion. - A free point particle of unit mass has the action 



Ao\x] 



/3 



dr i^(T). 



The amplitude (0/3|0 0)o is given by the Gaussian path integral 

(0/3|00)o= / Px(T)e-'^°[^l =e-(V2)^i°s(-^') = [27rp]-'^\ 



A coordinate transformation x{t) = /('z(t)) brings the action ( pO| ) to the form 



1 r 



dTg{q{T))q (r) 



where g{q) = f'^{q)- The measure 'Dx{t) = Yir^-'^i''') transforms as follows: 
Vx{t) = Y[dx{T) ^ J Yldqir) = JVqir), 



(40) 



(41) 



(42) 



(43) 



where J is the Jacobian of the coordinate transformation 



J 



^(1/2)5(0) /^"rfrlogg(q(r)) 



Thus the transformed path integral ( p| ) takes precisely the form 
the exponent 



Aot[g] = dr 

JQ 



(44) 

with the total action in 



lgiqir))eir)-^S{0)logg{q{r)) 



(45) 



This is decomposed into a free part 



1 



drq^ir) 



and an interacting part 

f/3 



(46) 



(47) 
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The path integral @ is now formally defined by the perturbation expansion 

(0,/3|0,0) = J 2?g(r)e'^''M-^-M = J I?g(T)e--^«M (^1 - Ant + ^Al, - . . . j 

= (2^/3)-i/2e-(A„t>c+i(^?„.>c-...^ 
with the harmonic expectation values 

(...) = (27r/3)i/2 /l?q(r)(...)e--^«M, 



(48) 



(49) 



and their cumulants {Af^^)c = (-^fnt) ~ (Ant)^, • ■ • containing only connected diagrams. If 
our calculation procedure respects coordinate independence, all expansion terms must van- 
ish to yield the trivial exact results As an example we shall consider the coordinate 
transformation 



= /(?) = 9 - ^£9^ + ^e"?^ 



such that 



(50) 



(51) 



where e is a smallness parameter. Substituting Eq. ( |5l[ ) into Eq. (|47|) yields up to second 
order in e: 



Ant [q] 



dr 



q\T)-SiO) 



(52) 



We shall now calculate order by order in e the expansion terms contributing to the square 
bracket in the second line of Eq. (|4J 



Diagrams. - To first order in e, the square bracket in the second line of Eq. (^) receives 
a contribution from the expectation values of the linear terms in e of the interaction (BS^ 



{A 



lin in ex 
nt / 



dT{eq^{T)e{T)-5{0)eq^(T)). 



(53) 



Thus there exists only three diagrams, two originating from the kinetic term and one from 
the Jacobian action: 

e OO +2£ CD ~em Q. (54) 

To order e^, we need to calculate only connected diagrams contained in the term {Af^^}/2 
in (|4^), all disconnected ones being obtainable from the cumulant relation (Af^^^) ~ {Af^^jc + 
{Aint)^- We distinguish several contributions. 

First, there are two local three-loop diagrams and one two- loop local diagram coming from 
the kinetic term and the Jacobian of the interaction (E3), respectively: 



3 ^ + 12 ^ - S{0) 00 
We call a diagram local if it involves only equal-time Wick contractions 



(55) 
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The Jacobian part of the action I p2^ contributes further the nonlocal diagrams: 

^ 1 2 52 (0) O - 4 <5(0) [ + CO + 4 00 ] } . (56) 

The remaining diagrams come from the kinetic term of the interaction ( p2[ ) only. They are 
either of the three-bubble type, or of the watermelon type, each with all possible combinations 
of the three line types (^) : The sum of all three-bubbles diagrams is 

^ [4 000 + 2 (XO + 2 CXDO + 16 aX) + 16 Oa:) + le 000 + le 000 ] • (57) 

The watermelon-type diagrams contribute 

(58) 



2! 



Path Integral in Curved Space. - Before we start evaluating the above Feynman diagrams, 
we observe that the same diagrams appear if we define path integral in a higher-dimensional 
target space q\ The generalization of the formal expression ( ^ ) is obvious: we replace g{q) 
by gij{q) in the kinetic term, and by g{q) —>■ det{gij{q)) in the measure, where gij{q) is the 
metric induced by the coordinate transformation from cartesian to curvilinear coordinates. 
In a further step, we shall also consider gij{q) more generally as a metric in a curved space, 
which can be reached from a flat space only by a nonholonomic coordinate transformation 
Jl9| . It was shown in the textbook that under nonholonomic coordinate transformations, 
the measure of a time-sliced path integral transforms from the flat-space form Yin '^^n to 
rin dq^/g^exp{AtRn/6), which has the consequence that the amplitude satisfies a Schrodinger 
equation with the pure Laplace-Beltrami operator in the kinetic Hamiltonian, containing no 
extra i?-term. Here we shall see that a similar thing must happen for perturbatively defined 
path integrals, where the nonholonomic transformation must carry the flat-space measure 



Vx^Vq^ exp / dr R/s] . (59) 




The proof of this rather technical issue is relegated to a separate paper. 

For n-dimensional manifolds with a general metric gij(^ we make use of the coordinate 



invariance to be proved by the vanishing of the expansion (54)-"(|5q). This will allow us to bring 
the metric to the most convenient Riemann normal coordinates. Assuming n-dimensional 
manifold to be a homogeneous space, as in a standard nonlinear cr-model, we expand the 
metric and its determinant in the normal coordinates as follows 

9^M) - + £ \ R^k,Jk, q'^'q"' + ^Rk^JkJ Rks^k,l q'^^ q^- q^^ q^^ + ..., (60) 

9{q) = l-£\R^3 q'q' + ^ Ru + \ i?™," Rkmi'^^ q'q'q^q' + ■■■ ■ (61) 

In our conventions, the Riemann and Ricci tensors are Rijk — diTj^ — . . . , Rj^. = Rjjk , and 
the curvature R = Ri has the positive sign for a sphere. The expansions ( |60| ) and (^l|) have 
obviously a similar power structure in as the previous expansion 
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In normal coordinates, the interaction (E^) becomes, up to order e^: 



Ant [q] 



q'q' 



drj e-Rikji q^q^ + e^—RmjJ Rrtsl q"'q''q'' 

Rkurq^q^q^q'}, 



e^SiO)R.,q^q^ +e'^6iO)R,rr./ 



(62) 



with the same powers of q^ as in Eq. (|52|). The interactions (|5^) and (|62|) yield the same 
diagrams in the perturbation expansions in powers of e. In one dimension and with the trivial 
vertices in the interaction ( |5^ ) , the sum of all diagrams will be shown to vanish in the case of 
a flat space. In a curved space with the more complicated vertices proportional to Rijki and 
Rij, the same Feynman integrals will yield a nontrivial short-time amplitude. The explicit 
-Rijfci-dependence coming from the interaction vertices is easily identified in the diagrams: all 
bubbles in (p6|)-(57) yield results proportional to R^j, while the watermelon-like diagrams ( p8| ) 



carry a factor R^jj^i- In our previous work 1^-^, all integrals were calculated in d dimensions, 
taking the limit c? ^ 1 at the end. In this way we confirmed that the sum of all Feynman 
diagrams contributing to each order in e vanishes. It is easy to verify that the same results 
are found using the procedure developed above. 

From Coordinate Independence to DeWitt-Seeley Expansion. - With the same proce- 
dure we now calculate the first two terms in the short-time expansion of the time-evolution 
amplitude. The results will be compared with the similar expansion obtained from the 
operator expression for the amplitude e^-^ with the Laplace-Beltrami operator = 
g-^/^dig^/^g'^(q)dj first derived by DeWitt M (see also M): 



{q,P\q',0)^{q\e^^'/^\q')^ 



1 



Y.P''ak{q,q'), 



(63) 



fc=0 



with the expansion coefficients being for a homogeneous space 



ao{q, q') 
ai{q,q') 
a2{q,q') 



1 

12 



R^J^q'^q' 



1 

360^ 



R\hR 



k l-^imjn 



1 
288 



RuR-n 



-R- 



12 
1 

288' 



144 



-R Ri 



1 
360 



R^^R, 



kilj ' 



1 

360 



ryklm D 



1 



180 



RlRnj Aq^Aq^ 



R^'^^^Rj^ki ~ H}"^ Rj] 

720 ^ 720 ^ 



where Aq^ = (g — q'Y . For Ag* = this simplifies to 



(0,/3|0,0) 



/2^ 



1 



12 72 



:R- 



10 



(R'^^^R 



ijkl 



R'^Rrj) 



(64) 



(65) 



The derivation is sketched in Appendix B. 

For coordinate independence, the sum of the first-order diagrams (|5^) has to vanish. An- 
alytically, this amounts to the equation 



/ dT [A(t, r)-A- (r, r) + 2 A ^(r, r) - 5{0)A{t, t)] 
Jo 



0. 



(66) 



In the cZ-dimensional extension, the correlation function A' (r, r) at equal times is the limit 
d ^ 1 of the contracted correlation function p,A^(a;, x) which satisfies the d-dimensional field 
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equation (|T]). Thus we can use Eq. ( |To| ) to replace 'A'(r, t) by (5(0) — 1/(3. This removes the 
infinite factor (5(0) in Eq. (pq) coming from the measure. The remain(ier is calculated (iirectly: 



dr 



A{t,t) +2A2(r, r) 



This result is obtained without subtleties, since by Eqs. (^) and (|^) 



whose integrals yield 



A(T,T)=r-— , ■A2(r,r) = -- 



1 /-/3 /■/3 
- dTA{T,T)= dT-A'(T,T) 



1 A(t,t) 



/3 

12' 



(67) 



(68) 



(69) 



The same first-order diagrams ( p4|) appear in curved space, albeit in different combinations: 



1 /"^ 

--i? / dr [A(r,r)A-(r,r) - A2(T,r) -5(0)A(r,r) 
6 Jo 



which is evaluated, using the integrals (p39|), to 



I dT 



iA(r,r) + A2(r,r) 



24 



(70) 



(71) 



This has to be supplemented by a similar contribution coming from the nonholonomically 
transformed measure (M). Both terms together yield the first-order DeWitt-Seeley expansion 



(0/3|0 0) 



1 



/2^ 



1 + ^^ 



(72) 



in agreement with (pq). 

We now turn to the evaluation of the second-order diagrams. The sum of the local diagrams 
(|55| ) is given by 

^(55) - -^£2 J^^ [3^2 ^) ^ i2A(r, r) A ^(r, r) - 6(0) A^r, r)] . (73) 
Replacing A'(r, r) in Eq. (fz^) again by (5(0) — and taking into account the equality 



/ (iTA(T, r) 

^0 



A(t, r) -4A2(t,t) 



= 



following from Eq. (bq), we find only the divergent term 



/3 



3(5(0) / dTA'(T,T) 



10 



,5(0) 



(74) 



(75) 



The sum of all bubbles diagrams (p6D-(57) resembles a Russian doll, where the partial 
sums of different diagrams are embedded into each other. Therefore, we begin the calculation 
with the sum (pq) whose analytic form is 



^^^^^ ^yJo lo '^^ {2'5'(0)A2(r, r') 

-4 (5(0) [A(r, r) A ^ (r, r') + 4 A (r, r) A(t, t') A (r, r') + A^ (r, r') A (r, r)] } . (76) 
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Inserting Eq. ( p^ into the last equal-time term, we obtain 

^ y^''^''drdr'{-2<52(0)A2(r,r') 

-45(0) [A(r, r) 'A ^ (t, t') + 4 A (r, r) A(r, r') A (r, t') - A2(t, t') //?] } . (77) 

As we shall see below, the explicit evaluation of the integrals in this sum is not necessary. 
Just for completeness, we give the result: 

j:,«) = ^{-...,o)|-«(o)[g+4|i-i.|]} 

We now turn to the three-bubbles diagrams (p7|). Among these, there exist only three 
involving the correlation function pAj^(a;, x') A(t, t') for which Eq. ( p^ is not applicable: 
the second, fourth, and sixth diagram. The other three-bubble diagrams in ( p7| ) containing 
the generalization fj^A^{x,x) of the equal-time propagator A'(t, r) can be calculated using 
Eq. (0). 

Consider first a partial sum consisting of the first, third, and fifth three-bubble diagrams 
in the sum (^7|). This has the analytic form 

E(57) = T r rrfTdT'{4A(r,r)A2(r,r')A-(r',r') 

1,3,5 ^ -^0 Jo 

+ 2 A-(T,T)A2(r,r')A-(r',r') + 16 A (t, t)A(t, t') A(r,r') A-(r',T')} . (79) 

Replacing A'(t, r) and A'(r',r') by 6(0) — 1/(3 we see that Eq. ( [79| ) contains, with opposite 
sign, precisely the previous sum ( [76| ) of all one-and two-bubble diagrams. Together the give 

f/3 rP 



5^(56) -f 5^(13) - '4 [ [ drdr'(~^Air,ryA'{r,T') 

1,3,5 ^ Jo Jo LP 



9 1 fi 

-t- ^ A2(r, r') - - A (r, r) A(r, r') A (r, r') } .(80) 



and can be evaluated directly to 



1,3,5 



By the same direct calculation, the Feynman integral in the seventh three-bubble diagram in 
(g) yields 

I7 = ^''^%rdr'A(r,T)A(r,T')A(T,r')A(T',T')=-^. (82) 



The explicit results (|81|) and (^2|) are again not needed, since the last term in Eq. (gO|) is 
equal, with opposite sign, to the partial sum of the sixth and seventh three-bubble diagrams 
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in Eq. (|57|). To see this, consider the Feynman integral associated with the sixth three-bubble 
diagram in Eq. (^: 

/e - r rdrdr'-A(r,r)A(r,r')-A-(r,r')A-(r',T'), (83) 
Jo Jo 

whose d-dimensional extension is 

= / d'^xd'^x' f,A{x,x)A{x,x'),,A,{x,x')A,{x\x'). (84) 

Jo "'0 



Adding this to the seventh Feynman integral ( p2[ ) and performing a partial integration, we 
find in one dimension 

p2 2 1-13 -If. 

5^(57) = y 16 (/6 + /7) - YJ, J, dTdT'jA{T,r)A{T,r')A{ry) 

^ , (85) 



2 V45 



where we have used c? [A (r, t)] /dr = —1/13 obtained by differentiating (68). Comparing 
(|85| ) with (80), we find the sum of all bubbles diagrams, except for the second and fourth 
three-bubble diagrams in Eq. (^7|), to be given by 

E(56) + E'(il) = y(-f^)- (86) 

The prime on the sum denotes the exclusion of the diagrams indicated by subscripts. The 
correlation function A' (r, r') in the two remaining diagrams of Eq. (^), whose d-dimensional 
extension is i^Ay(x,x'), cannot be replaced via Eq. (|l0|), and the expression can only be 
simplified by applying partial integration to the fourth diagram in Eq. (^), yielding 

rP I'll 

/4 = / / drdr' A(r,r)A(r,T')A-(r,r') A(r',T') 
Jo Jo 

d'^x d'^x' A(x, x) p,A{x, x')fj_A^{x, x')A^{x' , x) 



/ d'^x d'^x' A(x, x)A^{x', x')dl [^A{x, x')Y 
Jo 



Uo lo drdr'A{r,ryA{r',T')^[A\r,T')] 

The second diagram in the sum ( |57| ) diverges linearly. As before, we add and subtract the 
divergence 

rl3 nP 



h = I I drdr' A(r,r)A-2(T,T')A(r',r') 
Jo Jo 

rP rP 

= / / drdr' A{t,t)[-A^{t,t')~S^{t~t')]A{t',t') 
Jo Jo 

rP rP 

+ / / drdr' A''{t,t)S^{t-t'). (88) 
Jo Jo 
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In the first, finite term we go to d dimensions and replace (5(t — r') — > S{x — x') — —A,yi,{x, x') 
using tlie field equation (31). After this, we apply partial integration and find 



13 r0 



Jo 

13 rf3 



d'^xd'^x' A{x,x) l,Al{x,x')- Al^{x,x')] A{x',x') 

d'^x d'^x' {~^^, [A(x, x)] A^{x, x') f,A^{x, x')A{x',x') 

+ A{x, x)A,4x, x')Axx{x, x')dl [A{x' , x')]} 

rl3 rP 

/ / dT(ir'2{-A(T,r) A-(T,T')'A-(r,T')A(r',T') + 
Jo Jo 

A(r,T)A(T,r')A(r',T')A-(r,T')}. 



Jo 



(89) 



In going to the last fine we have used d[A(T,T)]/dT = 2 A (r, r) following from (|68|). By 
interchanging the order of integration r t', the first term in Eq. (p9) reduced to the 
integral (|87|). In the last term we replace A'(t,t') using the field equation (|) and the trivial 
equation 

dTe(T)(5(r) = 0. (90) 



Thus we obtain 
with 



T — _L rdiv 

h - U + ^2 



jR 
^2 



rdiv 

io — 



90 120 



7^ 
90 



/3 rl3 



Jo 



dTdT'A^{T,T)5^{T^T'). 



(91) 

(92) 
(93) 



Using Eqs. ( p7| ) and (|9l|) yields the sum of the second and fourth three-bubble diagrams in 
Eq. (|§): 



^(57) = ^(2/2 + 16/4) = £N 

2,4 I-'O "'0 



P rl3 

dTdT'A^{T,T)5^{T-T')+'^ 



(94) 



(95) 



Finally, inserting this into Eq. (86), we have the sum of all bubbles diagrams 
5^(57) + Y.^) - [ drdr'A^T, r)5\T " r') + ^| • 

Note that the finite part of this is independent of ambiguous integrals of type (p6| ) . 

The contributions of the watermelon diagrams ( p8|) correspond to the Feynman integrals 

rt3 rP 

^(58) 2e2 / / drdr' [A2(t, r') A-2(r, r') 
Jo Jo 

+ 4A(r,r')A(T,T')A(r,r')A-(r,r') + A2(t,t')A2(t,t')] • (96) 



The third integral is unique and can be calculated directly: 



rP rP 

/lo-/ dT dr' A2(r,r')A2(r,r') =£■ 
Jo Jo 



2^ 

90' 



(97) 
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The second integral reads in d dimensions 

d'^xd'^x A{x, x)f_iA{x, x')A,y{x, x) f_iAu{x, x). 
This is integrated partially to yields, in one dimension, 

/9 = -i/io-i j j dTdT'A{ry)A\TyyA{ry). 
The integral on the right-hand side is the one-dimensional version of 

Ig, ^ / d'^xd''x'A{x,x')Al{x,x') ^f,A{x,x'). 



(98) 



(99) 



(100) 



Using the field equation (|3f|), going back to one dimension, and inserting A(t, r'), A'(t, r'), 
and "A (t, r') from (^), (^), and (|^), we perform all unique integrals and obtain 



Inserting this and ( p7|) into Eq. ( |99| ) gives, finally, 

.1 



f 

48 



dr e^T)S{T) 



720 J 



(101) 



(102) 



We now evaluate the first integral in Eq. (96). Adding and subtracting the linear divergence 
yields 

Is= f r drdr' A'{T,T'yA'iT,T') 

JQ JQ 

= / drdr'A2(T,T') [A-2(T,T')"<52(T-r')] +£M / drdr' A^iT,T)S\T-TilD3) 
Jo Jo Jo Jo 

The finite part of the integral ( |l03| ) has the d-dimensional extension 

d''xd''x'A^{x,x') [^Al{x,x')-Al^{x,x')] (104) 
which after partial integration and going back to one dimension reduces to a combination of 



jR _ 



integrals Eqs. (h03) and (fOl): 



-2/9 + 2/9, 



dTe\T)S{T) + — \0- 



72 



The divergent part of /g coincides with /j'^ in Eq. (|93|): 



7-div 



P rP 



Jo 



drdr' A^{t,t)5'^{t ~ t') = 



(105) 



(106) 



Inserting this together with ( |97| ) and (102) into Eq. (|9^), we obtain the sum of watermelon 
diagrams 



^(58) = 2e^{h+Ah + ho) 

{ rP rP q2 fp 

= 2 dTdT'A'{T,r)S\T - r') - '-^ | 



dr e^{T)S{T) 



60 



(107) 
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For a flat space in curvilinear coordinates, the sum of the flrst-order diagrams vanish. To 
second order, the requirement of coordinate independence implies the vanishing the sum of 
all connected diagrams (|55|)-(p8). Setting the sum of Eqs. ([zsl), (|9^), and (107) to zero leads 
directly to the integration rule (2^) and, in addition, to the rule 



drdr' A^{t,t)6^{t-t')^5{0) / dTA^{T,T), 



(108) 



which we postulated before in Eq. ( pOD to cancel the S{0)s coming from the measure at the 
one- loop level. 

The procedure can easily be continued to higher-loop diagrams to define integrals over 
higher singular products of e- and i5-functions. In this way we obtain the confirmation of the 
rule (pO|). We have seen that at the one-loop level, the cancellation of S{0)s requires 



dTA{T,T)S{0) =5(0) / dTA(r,r) 



(109) 



The second-order equation ( |108| ) contains the second power of A(t, r). To n-order we find the 
equation 



J dn . . .dTnA{TuT2)d{n,T2) ■ ■ ■ A{Tn,n)d{Tn,n) = d{0) J drA^{T,T) 



which reduces to 



dTidr„ A"(ri,ri)52(^^ _ ^ ^(q) j dTA"{T,T), 



(110) 



(111) 



and this is satisfied given the rule (^) . See Appendix C for a general derivation of these rules. 

Let us now see what the above integrals imply for the perturbation expansion of the short- 
time amplitude in curved space in Riemann normal coordinates. Taking into account the 
nonzeroth contribution ( [72| ) of the first-order diagrams reproduces immediately the first term 
in the second-order operator expansion (p4): 



B_ 
12 



n2 

288^ 



(112) 



The sum of the local diagrams ( p5| ) involves both tensors R^.^ and R^jf^i- To order e^^ -^^e find 



30 V 36 



1 



-.Rt 



24 



:R 



'ijkl 



S{0) 
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1 



:Rt 



30 



tR 



ijkl 



(113) 



The contribution of all bubbles diagrams (pq) and (|57|) contains only R 



E(56) + E(i7l) = 



1080 ^"'^"'"' ^ 432 



2 /^^ d2 



■Km~e'^R 



(114) 



This compensates exactly the (5(0)-term proportional to R^j in Eq . (|113| ) and yields correctly 
the third second-order term —R^j/720 in the operator expansion (|64|). 

Before turning to the contribution of the second-order watermelon diagrams (58) which 
contain initially ambiguous Feynman integrals we make an important observation. Compari- 
son with Eq. (|64[) shows that Eq. (IIS) contains already the correct part of the second-order 
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DeWitt-Seeley coefiicient R^j^i/72Q. Therefore, the only role of contributions of the wa- 
termelon diagrams ( p8| ) which are proportional to must be to cancel a corresponding 
divergent part of the sum (113). In fact, the sum of the second-order watermelon diagrams 
( pq ) reads now, 

24 

!)• 



(/8-2/9+/10}, 



(115) 



where the integrals Is.h, and Iiq were given before in Eqs. (106), (105), (102), and 
Substituting these into Eq. (115) and using the rules (|2^) and (|^ we obtain 



E 



(58) 



24 



RIh f f drdr'A^ir, t)5\t - r') = e^^Rl,, <5(0), 



(116) 



thus compensating the (5(0)-term proportional to R^jj^i in Eq. (113) and no finite contribution. 

For one-component target space as well as for n-component curved space in normal co- 
ordinates, our calculation procedure using only the essence of the d-dimensional extension 
together with the rules ( pO| ) and (^6|) yields unique results which guarantee the coordinate in- 
dependence of path integrals and agrees with the DeWitt-Seeley expansion of the short-time 
amplitude. The need for this agreement fixes the initially ambiguous integrals /§ and /g to 
satisfy the equations 

120' 



4/q + 



(117) 



2/9-1-/1, 



0, 



(118) 



as we can see from Eqs. (107) and (115). Since the integral /lo = /3^/90 is unique, we must 
have /g = — /3^/720 and /|^ — — /3^/72, and this is what our integration rules indeed gave us. 

The main role of the dimensional extension in this context is to forbid the application of 
Eq. (p^ ) to correlation functions 'A'(r, r'). This would have immediately fixed the finite part 
of the integral /s to the wrong value — /9^/18, leaving only the integral /g which would 

define the integral over distributions (26). In this way, however, we could only satisfy one of 
the equations (117) and (118), the other would always be violated. Thus, any regularization 
different from ours will ruin immediately coordinate independence. 

It must be noted that if we were to use arbitrary rather than Riemann normal coordinates, 
one can fix ambiguous integrals already at the two-loop level, and obtains the conditions ( p^ ) . 
Thus, although the calculation in normal coordinates are simpler and can be carried more 
easily to higher orders, the perturbation in arbitrary coordinates help to fix more ambiguous 
integrals. 

Let us finally compare our procedure with the previous discussion of the same problem 
by F. Bastianelli, P. van Nieuwenhuizen, and others in Refs. [f7|-p5[. Those authors suggested 
for almost ten years two regularization schemes for perturbative calculation on a finite-time 
interval: mode regularization (MR) p|-pT|] and time discretization (TS) |11 13|. They gave 
a detailed comparison of both schemes up to three loops in Ref . . Their main goal was 
to calculate of trace anomalies of quantum field theory by means of path integrals P, [l^,[l5[]. 
From the present point of view of extended distribution theory, mode regularization (MR) 
amounts to setting 

dTe^{T)5{T) = \. (119) 



With this rule, the ambiguous integrals /§ and /g yield = — /3^/18, /g = /3^/180. However, 
these values do not allow for coordinate independence, nor do they lead to the correct short- 
time DeWitt-Seeley expansion of the amplitudes. This is what forced the authors to add 
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an unpleasant noncovariant "correction term" yl^^'^g'^ — — J dTT'^^^T\^^giig""-g^'^ /2A to the 
classical action, in violation of Feynman's construction rules for path integrals. In doing this 
they followed earlier work by Salomonson in Ref . |^ . 

Their time discretization scheme (TS), on the other hand, amounts to setting 

dTe^{T)5{T) =Q. (120) 

They applied this to purely one-dimensional calculations which, as we have shown in this 
paper, leads to the contradictory results depending on where partial integration or field equa- 
tions are used. While /§ is again = — /3^/18, the result for /g = 7/3^/360 is not unique. To 
obtain coordinate independence as well as the correct DeWitt-Seeley expansion, they had now 
to add another noncovariant "correction term" — J drV^^^Xlig^^ thereby following 

the original work of Gervais and Jevicki in Ref. . 

In recent papers [p^-|l8|, the authors of Refs. and [|l^ have begun following our 
method of dimensional regularization j2[ ^ , adapting it to a finite time interval in fl^ . They 
now obtain, of course, correct coordinate-independent results without noncovariant additional 
terms in the action. They do not, however, exhibit the precise location of ambiguities as we 
did here, and most importantly, they do not derive from their results rules for integrating 
products of e- and (5-functions, which are central to the present paper. In particular, they do 



not realize that dimensional regularization amounts to the integration rule Eq. (120) 
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Appendix A: Integrals I14 and //| from two-loop expansion in arbitrary coordinates. 
To order e, the metric and its determinant have the expansions: 



9i3 (9) 
logg(q) 

The interaction (EtI) becomes 



5ij + V£idkgtj)q'' + e^{^l^kg^j)q'q'', 



V^9''i^kg^J)q^ + e- g^= [{didm,) - 3""(5ig.™)(9fc5,n)]g'<z'. 



(121) 



AntM = / drj ]^^/e{dkgij)q^ +'^e{^l^kg^j)q^q^ 



q'q' 



1^6{0)g^^{dkg^,)q'' - y6{0)g^^[{dAg.,,) 



.g"" {digirn){dkgjr. 



q'q' 



}.(122) 

To the first-order in e, the perturbation expansion ( [48D with the interaction (122) consist 
of two sets of diagrams proportional to r,j^ k and F^^ j., respectively. First, there are the same 
l ocal diagrams as in Eq. (^^: the first two local diagrams coming from the kinetic p art o f 
(122) carry a factor Tij, k, while the last local diagram, coming from the measure part of (122), 
involves both factors Tij^ k and T^j f. . Omitting the T^j fc-part of the last diagram, the terms 
linear in the Christoffel symbol F^^ k coming from the sum of local diagrams in ( ^4|) reads 

^(54) = -|(a,9fcgy)^''dr{g^^5'='-A-(r,r)A(r,r)-f 25V' A ^(r, r)-5(0)g^^"/' A(t, r)} 
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(123) 



In addition, the interaction (122) generates nonlocal first-order diagrams proportional to 
Together with nonlinear in Christoffel symbol part of the last local diagram in Eq. (^^ , they 
are represented as follows 



■ e^i^ (5'' Ffc/" + 5^'= r.feO [ O-O - ^(0) ] 
- J(/ 5'=" r,,^" r,„,, + r,fc'= r,,' + 25*^" r,/ r.^') Q-O 



+ (r// r.fc-'" + r„fc^- r,,,,) [ >(-h 



-5(0) o: 



(124) 



The Feynman integrals associated with the diagrams in the first and second lines of Eq. (124) 
read 



111 
and 



dT dr' { A (r, r)A(r, r')A(T', r') - 2(5(0) A(t, t) A(t, r') + S^{Q) A(r, t')} (125) 



/i2= J jdrdr' {A(T,T)A(T,T')A(r',T') -,5(0) A(r,r)A(T,r')}: 



(126) 



respectively. Replacing in Eqs. (125) and (126) A(t, t) and A(T',r') by ,5(0) — 1//3 leads to 
cancellation of the infinite factors (5(0) and 5'^{0) coming from the measure, such that we are 
left with 



and 



/i2 = dr dr' A (r, r) A (r, r') = - 



12' 



(127) 



(128) 



The Feynman integral of the diagram in the third line of Eq. (124) has (i-dimensional 
extension 



/l3 = 



dr dr' A (t, t) A (r' , t') A (t, t') 



(129) 



Integrating this partially yields 



/i3 = ^ y j dT dr' A(t, r') a (t', t') ^ ^ j^^ dr' A (r, r) A (r, r' 



12' 



(130) 



where we have interchanged the order of integration r ^ t' in the second line of Eq. (130) 
and used d[A(T, T)]/dr = -1//3. Multiplying the integrals ( p7| ), ( [l28| ), and ( |T30| ) by corre- 
sponding vertices in Eq. (124) and adding them together, we obtain 



^(i24)^|^.9V'r./r, 



1,2,3 



24 



(131) 
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The contributions of the last three diagrams in the fourth and the fifth fine of Eq. (124) 
correspond to the ambiguous integrals (p ) and jl^), respectively. Moreover, the difference 
of two diagrams in the last line of Eq. (124) contains only the finite part of the integral ([T^, 
since its divergent part (|^) is canceled by the contribution of the local diagram with the 
factor (5(0). Multiplying these integrals by corresponding vertices in Eq. (124) yields the sum 
of diagrams in the fourth and the fifth line of Eq. (124) as follows 

^(124) = I [g'"^ g^' E,," T^j, „ (/u + /fs) + 5'' ^ij' ^^k' (3/i4 + I^) }■ (132) 

4,5 



On the other hand, to guarantee the coordinate independence of path integrals, this sum 
must be 



^(i24) = -|f5Vr.,"r,,.„. 



4,5 



24 



(133) 



Adding this to (131), we find the sum of all diagrams in (124) as follows 



^(124) - ^5^^"/' (r./r,i,„ - E,fc"E,7.„) 



24 



(134) 



Together with the sum over all diagrams in (54) calculated in (123) this yields, finally, the 
sum of all first-order diagrams 



(135) 



The result is perfectly covariant and agrees, of course, with Eq. ( [7l|) derived in normal coor- 
dinate. Comparing now Eq. (132) with (|l33|), we find 



3/1 



rU 

^15 

tR. 
-'15 



12 ' 



(136) 



Thus, coordinate independence specifics the initially ambiguous integrals (|l]) and (|l|) to 
have indeed the values (p^). 

Appendix B: Operator derivation of short-time DeWitt-Seeley expansion. - Here we give 
a short derivation of the DeWitt-Seeley expansion ( |63| ) . In a neighborhood of some arbitrary 
point (7q we expand the Laplace-Beltrami operator in normal coordinate system ( |6l| ) as 

D^ = d^- ^-R^k^jk, {qo){q - qo)"' {q - qof'd.d, - '^-R^,{q^){q - qofdj. (137) 

To find the coefficients ak{q,q') in Eq. (|6^), we resort to perturbation theory. The time 
displacement operator H = in the exponent of Eq. (^3|) is separated into a free part 

Hq and an interaction part Hmt as follows 



Ha 

Hint 



-^Rik^jk^ {q - qof^ {q - qof^^^^j 4- -Rij{q - qofdj 



(138) 
(139) 
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The transition amplitude (p3) satisfies the integral equation 



(g,/3 I q',0) = {q I e-^(^°+^'"') | q') = {q \ e"^^" 

f/5 



where 



1 - / dae""" Hinte 
Jo 

= {q,P\q', 0)o - £da J d^q [q, P - a \ q, 0)o i?int (g) [q, a | g, 0) 

1 



(g,/3|g',0)o-(g|e-^«° U') - 



,-(A9)V2/3 



To first order in 7?int we obtain 

(q, P I g', 0) = (g, /3 I g', 0)o - J^a J d"q (g, /3 - a | g, 0)o i?i„t (g) (g" a | g, 0)o. 
Inserting (|139D and choosing qo = q', wc find 



(140) 
(141) 

(142) 



(g,/3|g',0) = (g,/3|g',0)oa + 



1 







1" 









d"(Ag) 



-[Ag-(<T/,3)Ag]V2a 



--R.k.jk.Aq'^Aq''- ( ~ — 



(5^^ Ag^AgA 1^ Ag^Ag 



(143) 



where we have replaced the integrating variable g by Ag = g — g' and used the notation 
a = (/3 — a)a/p. There is initially also a term of fourth order in Ag which vanishes, however, 
because of the antisymmetry of Rikji in ik and jl. The remaining Gaussian integrals are 
performed after shifting Ag Ag + ct Ag//3, and we obtain 



(g,/3|g',0) = (g,/3|g',0)o{l + -y^ da[^i?,,(g')Ag*Ag-' + ^i?(g')] } 
= (g, /3 I g', 0)o [l + (g') Ag^ Ag^' + ^R{q') 



(144) 



Note that all geometrical quantities are evaluated at the initial point g'. They can be re- 
expressed in power series around the final position g using the fact that in normal coordinates 



5y(g') = .gy(g) + ^R.k,jkAq)Aq'''Aq''-' + 

5.,(g')Ag*Ag^ -g.,(g)Ag^Ag^ 



(145) 
(146) 



the latter equation being true to all orders in Ag due to the antisymmetry of the tensors Rijki 
in all terms of the expansion (145), which is just another form of writing the expansion (pG) 
up to the second order in Ag*. 

Going back to the general coordinates, we obtain all coefficients of the expansion ( |63| ) 
linear in the curvature tensor 



(g,/3|g',0) 



-g,j(9)Ag'AgV2/3 



1 



(3 



l + -:R^AQWAq^ + ^Riq) 



12 



12 



(147) 



The higher terms in d63) can be derived similarly, although with much more effort. 
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A simple cross check of the expansion ( |63| ) to high orders is possible if we restrict the space 
to a sphere of radius r in D dimensions. Then 



R 



1 



ijkl 



i9ik 9ji- gii9]k), i, j = 1, 2, . . . , n = £> - 1, 



(148) 



where n = 2 is dimension of a sphere, and D = 3 is dimension of a flat embedding space, 
respectively. Contractions yield Ricci tensor and scalar curvature 



Ri 



Rki 



D -2 



R = Ri^ — 



{D-l){D-2) 



and further: 



r2 _ 



2 (D -!)(£)- 2) 



, {D-l){D-2f 



(149) 



(150) 



Inserting these into (Q) , we obtain the DeWitt-Seeley short-time expansion of the amplitude 
from g = to g = up to order 0^: 



(0,/9|0,0)- 



/2^ 



D-l 



1 + {D-l){D-2)±-^ + iD-l){D-2){5D' 17D + 18)^ 



(151) 



On the other hand, we may follow Ref. p2| , and calculate explicitly the partition function 
for this system 



Z{P) = ^ exp[-l{l + D-2)/3/2r2] 



(152) 



1=0 



where + D — 2) are the eigenvalues of the Laplace-Beltrami operator on a sphere and 
di = (2^ + £' — 2)(Z-|-£' — 3)!/Z!(r' — 2)! their degeneracies. Since the space is homogeneous, the 
amplitude {0,(3 \ 0,0) is obtained from this by dividing out the constant surface of a sphere: 



(0,/3|0,0) 



(153) 



For any given D, the sum in (152) easily be expanded in powers of (3. As an example, take 
D = 3 where 



Z{(3) = ^ (2/ + 1) exp[-l{l + l)/3/2r2] 



(154) 



1=0 



In the small-/? limit, the sum (154) is evaluated as follows 

/•oo 

Z{I3) = d + 1)] exp[-?(/ + l)/3/2r2] + ^ (2^ + 1) [l - l{l + l)p/2r^ + . . . ] . (155) 

1=0 



The integral is immediately done and yields 



dz exp(— z/3/2r^ 



2r2 



(156) 
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The sums are divergent but can be evaluated by analytic continuation from negative powers of 
I to positive ones with the help of Rieniann zeta functions (^{z) — X^i^i ' 
for all even negative arguments. Thus we find 

oo oo 

^(2^ + 1) = l + ^(2/ + l) = l + 2C(-l)-- = -, 



1=1 



2r2 



^(2Z + l)/(/ + l) 



1=0 



2/' 



1=1 



2r2 



[2C(-3) + C(-l)] 



which vanishes 

(157) 
(158) 



Substituting these into (155), we find 

2r2 



6r2 60r4 



(159) 



Dividing out the constant surface of a sphere Attt^ as required by Eq. ( |153[ ) , we obtain indeed 
the expansion (151) for D = 3. 

Appendix C: Cancellation of all powers of 6(0). - There is a simple way of proving the 
cancellation of all UV-divergences S{0). Consider a free particle whose mass depends on the 
time with an action 



Aot [q] 



dr 



iz(r)g2(r)- 1^(0) log Z(r) 



(160) 



where Z{t) is some function of r but independent now of the path (j(t). The last term is 
the simplest nontrivial form of the Jacobian action in (^). Since it is independent of g, it is 
conveniently taken out of the path integral as a factor 



J 



^(1/2)5(0) /^"dr log Z(r) 



We split the action into a sum of a free and an interacting part 

Ao= J^dT^q'ir), Ant- dT^[Z{T)-l]eiT), 



(161) 



(162) 



and calculate the transition amplitude (48) as a sum of all connected diagrams in the cumulant 
expansion 



(0,/5|0,0) = J y I?(z(T)e^-^°M--^""[''l = J y I?g(r)e--^«M 



1 - Air 



Ia - 

2 



(27r^)-i/V 



int/ 



liAl 



(27r/3)-i/Ve-<-^'"'>^+^<-^-'> = 



(163) 



We now show that the infinite series the of (5(0)-powers appearing in a Taylor expansion of 
the exponential (161) is precisely compensated by the sum of all terms in the perturbation 
expansion (162). Being interested only in these singular terms, we may extend the r-interval 
to the entire time axis. Then Eq. ( p^ yields the propagator 'A'(t, t') = (5(t — r'), and we find 
the first-order expansion term 



(Ant)c = J dT^[Z{T)-l]-/^{T,T) = -^S{0) J rfr[l-Z(r)]. 



(164) 
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To second order, divergent integrals appear involving products of distributions, thus requiring 
an intermediate extension to d dimensions as follows 

(A'„t>c - J I dndT2^{Z -1),^{Z -l)22-A{Tur,)-A{T2,n) 

d'^Xi d'^X2 ^{Z - 1)1 i(Z - 1)2 2 f,A^{xi,X2) ^A^,{X2,X1) 

d'^xi d''x2 \{Z - 1)1 i(Z - 1)2 2 A^^(x2, xi) A,,[x^,x2) , (165) 

the last line following from partial integrations. For brevity, we have abbreviated [1 — Z(Ti)\ 
by {1 — Z)i. Using the field equation ( |3l| ) and going back to one dimension yields 



To third order we calculate 



dTldT2il-Z)l{l-Z)2S^{Tl,T2). 



(166) 



{A 



int/c - 



dT2 dT3 i (Z- 1) 1 i (Z- 1)2 i (Z- 1)3 8 (n , T2 ) A {t2 , T3) A (rg , ri ) 

d'^xi d'^X2 d'^xs \^{Z-\)i \^{Z-\)2 \{Z-\)-i 8 ^A^(a;i, 0:2) ^A<,(a;2, xg) ^A^(x3, xi) 

d'^xxd'^X2d'^x:!,h^Z-\)^ 1)2 ^(Z- 1)3 8 A^^(x3, xi) A^^(a;i, X2) A„„(x2, xg). 

(167) 

Applying again the field equation ( |3l| ) and going back to one dimension, this reduces to 



(•^fiit)c 



dri dT2 dT3 (1 - Z)i (1 - Z)2 (1 - Z)3S{ti,T2) S{t2,T3) 5(t3, ri). (168) 



Continuing to n-order and substituting Eqs. (164), ( |166| ), (168), etc. into (163), we obtain in 
the exponent of Eq. (163) as sum 



(^int)c + r, (•^int)c gj (-^int/c ' 



00 



n 



with 



where 



dTi . . . dT„ C{ti,T2) C{t2,T^) . . . C{Tn, Ti) 



C{t,t') = [1-Z{t)]6{t,t'). 



Substituting this into Eq. (170) and using the rule ( pO| ) yields 

^n = J J dndTr. [1 - Z{n)r6^{n - t„) = S{0) j dT[i- Z{t)\^. 

Inserting these numbers into the expansion ( |169| ), we obtain 



(169) 

(170) 
(171) 

(172) 



(A, 



--5(0) I dr log Z{t), 



(173) 



which compensates precisely the Jacobian factor J in (|163|) 
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